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We investigate the properties of dibaryons containing u and d quarks in the constituent quark 
model. In constructing the ground state wave function, we choose the spatial part to be fully 
symmetric and the remaining color, isospin and spin part to be antisymmetric so as to satisfy the 
the Pauli principle. By adapting the IS coupling scheme that combine the isospin basis function 
with the spin basis function, and subsequently coupling this to the color singlet basis function, we 
construct the color ® isospin ® spin states compatible with the physical states of the dibaryon. 

By using the variational method, we then calculate the mass of the dibaryon in a nonrelativistic 
potential model, involving Coulomb, color confinement and color-spin hyperfine interaction. In 
particular, to asses the stability for different types of the confinement potential, we introduce one 
that is linearly proportional to the interquark distance and another to its square root. For all cases 
considered, we find that there are no compact bound states against the strong decay. 

PACS numbers: 14.40.Rt,24.10.Pa,25.75.Dw 


I. INTRODUCTION 

The recent obervation of d*(2380) with quantum num¬ 
ber J(J P )=0(3 + ) mesasured by the WASA detector at 
cosy G® revived interests in the study of multiquark 
hadrons, and led to the renwed investigation of the pos¬ 
sible existence of either AA or six-quark state. Theo¬ 
retically, starting from the work of Jaffe 0, there were 
already many studies on the stability of a six-quark sys¬ 
tem. In particular, in relation to d*(2380), a study of 
nonstrange diybaryon was made in Ref. and a work 
in Ref. [9| was based on one gluon-exchange interac¬ 
tion. Initially, using the Bag model with strange quark, 
Jaffe (?) predicted that the H-dibaryon with J p =0 + and 
1=0 consisiting of uuddss could be stable against de¬ 
cay into two A baryons, when only the color-spin hyper¬ 
fine interaction was taken into account. Using a similar 
quark model, Silvetre-Brac and Leandri m classified all 
dibaryon states within the SU(3)f representation, and 
investigated the stability of these states against the de¬ 
cay into the allowed two baryon decay: Through this 
study, they found that the UU dibaryon is most likely 
bound, and H-dibaryon could be stable. 

Additional models were used to study the stability 
of the H-dibaryon; these include lattice gauge a bag 
model IT3 . Skyrme model ITdl ] , and potential model Ga¬ 
in Ref. [l5j], using the Goldstone boson exchange inter¬ 
action, the authors predicted that H-dibaryon could not 
exist. While experimental searches for H dibaryon seems 
to suggest that it is not stable against strong decay for 
realistic quark masses, recent lattice gauge theory cal¬ 
culations suggest that it does become bound when the 
quark mass increases fltil fl7j . 

In addition to the study of H-dibaryon, the dibaryon 
with strangeness -1 or -3 has been proposed by Malt- 
man (is! a nd Goldman |]j§], respectively. Pepin and 
Stancu [20] investigated the stability of uuddsQ (Q=c 
or b ) type of quark configuration with a chiral con¬ 
stituent quark model, and found the dibaryon to be un¬ 


stable against strong decay. The stability of a multiquark 
system is known to increase when heavier quarks are in¬ 
cluded in the tetraquark configuration ( qqbb ) fjll [2j| . Al¬ 
though the mechanism for stability is different, dibaryons 
with heavy quarks, such as q A Q 2 (Q=c or b), have been 
studied within the simple chromomagnetic model [25| . 

In this paper, we investigate all dibaryon states con¬ 
taining the u and d quarks, and calculate their masses in 
the frame work of nonrelativistic quark model by using 
the variational method, with a potential that includes the 
color spin hyperfine potential introduced in Ref. [26| . In 
order to examine the stability of the dibaryon for strong 
decay, we first fit the parameters of the model to repro¬ 
duce the masses of the baryon multiplet. Then, by com¬ 
paring the dibaryon masses to the relevant two baryons 
threshold, we determine whether the dibaryon is bound 
against strong decay. 

The confinement part in our model originates from the 
effect of one gluon exchange interaction A^A^. But in 
principle the Hamiltonian in the SU(3) symmetric quark 
model can also have a term proportional to the SU(3) 
invariant operators in the cubic form which could origi¬ 
nate from an intrinsic three-body color confinement in¬ 
teraction. There are two independent three-body color 
invariant operators; one that can be expressed in terms 
of two different types of Casimir operator of SU(3), the 
other that can not. Since we expect that adding a three- 
body color invariant operators is very important for the 
stability of a dibaryon with different flavors, it is neces¬ 
sary to introduce a formula for the operators in terms 
of the element of the permutation grou p Sq , based on 
the established-formula by Stancu [27| and Dmitrasi- 
novic [Hi]. Using these formula, we calculate the matrix 
element of the operators with respect to the color singlet 
basis function in six quarks system, and explore the role 
of the operators to baryon and dibaryon masses. 

This paper is organized as followings. We introduce 
the Hamiltonian and fit the baryon spectrum in section 
II. We construct the spatial function in a simple Gaussian 


2 


form in section III. We present all of the physical states 
and construct the color (g> flavour ® states of dibaryon 
in section IV. We show the numerical results obtained 
from the variational method, and deal with the three- 
body color operators in section V. Finally, we give the 
summary in Section VI. 


II. HAMILTONIAN 


For the nonrelativistic Hamiltonian, we take the con¬ 
finement and hyperfme potential for the color and spin 
degree of freedom given by 


i= 1 i<j 


C 

ij 


+ V? s ), (1) 


where Wj’s are the quark masses, and A^/2 are the color 
operator of the Pth quark for the color SU(3), and V£j 
and V$ s are the confinement and hyperfme potential, 
respectively. For the confinement potential, we adopt 
the following two different types : 


• Type 1 ; 


V c = 

*7 




( 2 ) 


In the following analysis, we take the units for k, 
a 0 , and D to be MeVfm, (MeV) -1 fm, and MeV, 
respectively. 


• Type 2 ; 


V-- =- 

1"ij 


( r b ') 1/2 


- D. 


a 0 


(3) 


Here, the units of a o is taken to be (MeV) 1 (fm) 1 / 2 . 


The hyperfme term which effectively splits the multiplets 
of baryon with respect to spin is given by, 


V? s 

u 


fr 2cV 1 c -(nH 2 /(r w )V . fj. 

mirrijC 4 (r 0 ij) 2 rij 1 3 ' 


(4) 


where the unit of n' is taken to be MeVfm. Here, r,; ; - is 
the distance between interquarks, | r, — r j |, and (raj) 
are chosen to depend on the masses of interquarks, given 

by, 

i // | n \ / c \ 

roij = l/(« + P ---), (5) 

rrii + rrij 


where the unit of a is (fm) -1 , and the unit of /3, 
(MeVfm) -1 . We choose to keep the isospin symmetry 
by requiring that m u =rrid (MeV). In the Hamiltonian, 
the parameters have been chosen so that the fitted mass 
of both baryon octet and decuplet are comparable with 
those of experiments. 


TABLE I: Parameters fitted to the experimental baryon octet 
and decuplet masses for the two different types of potentials. 



Hi 

ao 

D 

k! 

a 

P 

m u 

m s 

Type 1 

107.6 

0.001062 

952.6 

107.6 

2.36 

0.0015 

340 

610 

Type 2 

109.6 

0.001103 

963.6 

168.6 

2.16 

0.0018 

348 

612 


A. Baryon Spectrum 


In constructing the basis function of a baryon, we re¬ 
strict the flavor symmetry to isospin part only, so that 
we consider only the u, d quarks as identical quarks, and 
find the total wave function of baryon, according to Pauli 
principle. When we calculate the expectation value of 
the potential terms for baryon with certain symmetry, 
it is convenient to introduce the following three Jacobian 
coordinates so as to reduce our problem to the two body- 
system in the center of mass frame : 

• Coordinate I ; 


Xl = 7f( ri_r2) ’ 


x 2 = 


(6) 


Coordinate II 


* 1 = 7§ (ri-rs> ’ 


2/2 = 


2(r 2 - I ri - irj). (7) 


Coordinate III 


* = Tfi ' 2 - r3> ’ 


Z2 = 


2 1 1 . , . 
j(ri- 5 r 2 --r 3 ). ( 8 ) 


By using simple Gaussian function, we construct the 
following fully symmetric spatial function for baryons 
composed of u and d constituent only : 

R = exp[-a(xi) 2 - b(x 2 ) 2 } + exp[-a(j/i) 2 - b(y 2 ) 2 }+ 
exp[—a(zi) 2 — b(z 2 ) 2 ], (9) 

where a and b are variational parameters. Since the total 
wave function of baryon, such as N and P (1=1/2, S=l/2) 
or A (1=3/2, S=3/2) is fully antisymmetric due to the 
Pauli principle, the rest of the total wave function must 
be fully antisymmetric, if we choose the spatial function 
to be fully symmetric. 

Concerning the color basis function of the baryon, we 
consider the color singlet state as the hadrons are ob¬ 
served to be colorless. The baryon has only one color 
singlet state, coming from the irreducible representation 
of [3]c ® [3]c ® [3]c, given by, 

= ^W(iy'(2)g fe (3)), (10) 


I C)= 
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which is fully antisymmetric under the exchange of any 
two particles among 1, 2, and 3. We note that the 
Young tableau follows the rule of the standard Young- 
Yamanouchi representation, which will be shown later in 
detail. 

For the spin basis functions, the bar 
S=l/2 consisting of two different types, ai 
taining one type as follows : 


S=3/2 ; 

|S 3 / 2 )= 


=ttt- 


(23) 


1 

2 

1 
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( 11 ) 


Now, we can construct the basis function of the color 
® isospin ® spin for (1=3/2,S=3/2) and (1=1/2,S=l/2) 
that are completely antisymmetric. These are given by 

I/ 3 / 2 , S 3/2 ) = 1C) 0 |T 3 / 2 , S' 3 / 2 ), 


• S=l/2 ; 

l^ 1/2 )= 

| 5 2 1/2 )= 


1 

2 

3 


1 

3 

2 



=75 (2 m-nt-m 


As we can see, the spin part of the basis functions of 
baryon IS 3 / 2 ) is completely symmetric, while that of 
| S\^ 2 ) and |S 2 /2 ) are partially symmetric; the former be¬ 
ing symmetric between particle 1 and 2, and the latter 
antisymmetric between particle 1 and 2. 

Likewise, we construct the isospin basis function of 
baryon for 1=3/2 and 1=1/2 : 


• 1=3/2 ; 
|/ 3 / 2 ) = 


1 2 


3 


=uuu. 


• 1 = 1/2 ; 

I TV 

iC>= 


1 

2 

3 


1 

3 

2 



= 7^(2 uud — udu — duu ), 
=-±=(udu — duu). 


|J 1/2 , S 1/2 ) = |C) ® |/ 1/2 , S 1/2 ). (12) 

For the hyperon, we treat the strange quark as distin¬ 
guishable from the u and d quarks and will not require 
the total wave function to be fully antisymmetric. It is 
then easy to construct the total wave function for the hy- 
perons with strangeness s=-l and s=-2, that satisfy the 
Pauli principle in the u and d quark sectors only. 


TABLE II: This table shows the mass of baryons in octet and 
decuplet obtained from the variational method. The fourth 
row indicates the experimental data. ( unit ; MeV ) 


(I,S) 

(I I) 

N, P 

(I 

v 2 5 2 / 

(0,|) 

A 

(b|) 

E 

(I 1) 

(1.1) 

E* 

(3 3, 

A 

Type 1 

977.1 

1315.3 

1115.6 

1206.0 

1530.2 

1403.4 

1267.5 

Type 2 

976.3 

1380.3 

1115.6 

1238.2 

1593.0 

1419.3 

1237.2 

Exp 

938.2 

1314.8 


1189.3 

1530 

1382.8 

1230 


939.5 

1321.7 

1115.6 

1197.4 

1531.8 

1387.2 

1234 


For the baryons with (1=3/2,S=3/2) and 
(1=1/2,S=l/2), the antisymmetry property of the 
total wave function can be easily obtained from the 
direct product of the totally antisymmetric part of 
the color basis function times the totally symmetric 
properties of the remain function comprised of the 
spatial function, isospin basis functions, and spin basis 
functions. The fully symmetric part of the isospin ® 
spin basis function in both cases of the isospin and spin 
and can be written as 

• A (1=3/2,S=3/2) ; 

|/3/ 2 , S , 3/ 2 ) = |/3/ 2 ^ g) l^ 3 / 2 ), 


. N, P (1=1/2,S=l/2) ; 

| / 1/2, 5 1/2) = _1 5( | Ji 1/ 2 ) g, | 5 V2 ) + | J l/2 ) g, |5 2 V2 )) . 

In the case of (1=1/2,S=l/2), we can see that 11 1 / 2 , S' 1 / 2 ) 
is symmetric between particle 1 and 2. The remaining 
symmetry for the permutation (23) can be deduced from 
the following formulas, according to the rule to the stan¬ 
dard Young-Yamanouchi representation : 


III. SPATIAL FUNCTION 


In order to construct a fully antisymmetric wave func¬ 
tion of dibaryon containing only identical u and d parti¬ 
cles, we choose the spatial function to be fully symmetric 
such that the rest of the wave function represented by 
color ® flavor ® spin should be antisymmetric. Since 
we restrict the SU(3) flavor of dibaryon to the isospin 
symmetry only, the flavor state can be identified with 
the isospin quantum number. In describing the system 
consisting of six quarks, it is convenient to deal with the 
system in the center of mass frame, reducing the number 
of suitable Jacobian coordinates of the system to five. 
The five Jacobian coordinates are given by 


®i = ^=(n-r 2 ), x\ = i(r 3 -r 4 + r 5 -r 6 ), 


x \ = ~ r4 ~ r 5 + r e), 


-(r 3 + r 4 -r 5 -r 6 ), 


x, = 


\/l2 


(r 3 + r 4 + r 5 + r 6 - 2r 4 - 2r 2 ). 


(13) 
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The variational method for calculating the mass of 
dibaryon turns out to be easy when the Gaussian form 
with respect to the Jacobian coordinates are used for the 
spatial wave function. Using the Jacobian coordinates 
given in Eq. (1131) in the Gaussian wave function, we find 
the form to be symmetric under the exchange of any two 
particles among 3, 4, 5, 6 , and at the same time symmet¬ 
ric under the exchange of two particles between 1 and 2 ; 
these symmetry properties are denoted as [3456] [12]. In¬ 
troducing the variational parameters a, b , c, the spatial 
function is then given by 

R Sl = exp[-(a(xi ) 2 + b(x \) 2 + b(x\) 2 + b(x\) 2 + 

c(^) 2 )]. (14) 


In addition to this Gaussian function, the full symme¬ 
try of the spatial function requires the linear sum of 14 
additional Gaussian functions, each of which has a spe¬ 
cific symmetry under particle exchange. The next set of 
five Jacobian coordinates with the symmetry of [2456] [13] 
is given by 


®i = ^(ri-r 3 ), x\ = i(r 2 - r 4 + r 5 - r 6 ), 

x l = ^( r 2 - r 4 - r 5 + r 6 ), x\ = i(r 2 + r 4 - r 5 - r 6 ), 


4 = 


\[V2 


(r 2 + r 4 + r 5 + r 6 - 2ri - 2r 3 ). 


(15) 


The corresponding Gaussian function specifying the sym¬ 
metry of [2456] [13] is given by 

R S2 = exp[—(a(x 2 ) 2 + b(x \) 2 + b(xl) 2 + b(xl) 2 + 

cOrl) 2 )]. (16) 


We find that the set of Jacobian coordinates necessary 
to obtain the fully symmetric wave function under the 
exchange of any two particles among 1, 2, 3, 4, 5, and 6 , 
and consequently the corresponding Gaussian functions, 
are the ones with the following symmetry; [3456] [12], 
[2456] [13], [2356] [14], [2346] [15], [2345] [16], [1456] [23], 
[1356] [24], [1346] [25], [1345] [26], [1256] [34], [1246] [35], 
[1245] [36], [1236] [45], [1235] [46], [1234] [56], 

Combining these Gaussian functions with the symme¬ 
try into a linear form, we obtain the spatial function with 
three variational parameters a, b , c which is fully sym¬ 
metric as follows 


R s =R S i + R s * + R s * + R s 4 + R°s + 

R Se + R s7 + R s8 + R s9 + R S10 + 

R S11 + R Sl2 + R S13 + R Sli + R S15 . (17) 

It is easy to check the symmetry of the spatial function 
with respect to all of the permutation of Sq, by consid¬ 
ering only the 5 permutations (12), (23), (34), (45), and 
(56), as these permutations generates all of the permuta¬ 
tions of Sq. 


IV. CLASSIFICATION OF DIBARYON WITH 
ISOSPIN SYMMETRY 


A. Isospin and spin state of dibaryon 


In this section, we investigate the state of the dibaryon 
consisting of identical u, d quarks, whose flavor part is 
characterized by isospin symmetry. Since the color ® 
isospin ® spin state of each quark can be represented 
by [3]c 0 [ 2 ]/ ® [ 2 ]g, the direct product of six quarks 
enable us to classify all of the states of dibaryon with re¬ 
spect to the state of isospin and spin, denoted by ]/, S). 
In our notation, the [3]c indicate the fundamental rep¬ 
resentation of SU(3)c, [ 2 ]/ the fundamental represen¬ 
tation of SU(2)j , and [2]g the fundamental representa¬ 
tion of SU(2)s■ In our case where we choose the spatial 
function of dibaryon to be fully symmetric, the color CS> 
isospin ® spin state of the dibaryon will be chosen to be 
fully antisymmetric. The fully antisymmetric state of the 
color ® isospin ® spin state can be easily obtained from 
the classifying of the multiplets of the direct six product 
of [ 12 ] C is, which is the fundamental representation of 
SU(12)cis, and gives the multiplet with dimension 924 
represented by Young tableau [l 6 ]. Using the original 
representation of [3] <7 ® [4]/s, which we will equivalently 
represent as ([3]c, [4]/s), the totally antisymmetric mul¬ 
tiplet of [l 6 ]c/S can be decomposed as 

[l 6 ]c/s =([l]c, [50]/s) © ([ 8 ]c, [64]/s) 0 ([10] c , [10] /s )© 
([1 O]c,[1O]/s)0([27]c,[6]/s). (18) 


By using the Young tableau, we can easily find that the 
multiplets in the right hand side of Eq. (18) is fully an¬ 
tisymmetric. According to the permutation group the¬ 
ory, for a given Young tableau, the fully antisymmetric 
function can be constructed by multiplying the Young 
tableau by its conjugate of the Young tableau, where the 
conjugate representation of a given Young tableau can 
be obtained by exchanging the row and column in the 
Young tableau : 


([!]c, [50]/s)= 



([ 8 ]c, [64] /s )= 
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me, [ 6 ]/ S )= 



Because the dibaryon is supposed to be a physically 
observable color singlet state, the dibaryon belongs to 
the five independent color singlet states represented by 
the Young tableau of [2,2,2] from the classification of the 
multiplets of the direct six times of [3]c- Hence, only 
([l] c , [50]/s) states in Eq. (1TH1) is allowed as the physical 
states of the dibaryon. We can also find the decompo¬ 
sition of [50] is with respect to the multiplets of [2]/ 0 
[2]s, given by, 

[50]/s =([1 ]/, [3]s) ® ([3]/, [l]s) ® ([3]/, [5]s)® 

([5]/i [3]s) © ([7]/, [l]s) © ([!]/> [7]s)- (19) 

These |J, S) states are all the possible states of dibaryon 
with isospin symmetry which satisfy the antisymmetry 
property and the color singlet requirement as an observ¬ 
able. As we see in the Eq. in, the isospin and spin 
of the dibaryon with isospin symmetry has I=0(S=0), 
I=1(S=1), I=2(S=2), and I=3(S=3), coming from the 
classification of SU(2) of the dibaryon, given by, 

[2]/ ® [2]/ ® [2]j 0 [2]/ 0 [2]/ 0 [2]/ = 

[l]i=0 0 f[3,3] © [3]/=l 0 -F[4,2]® 

[5] 1=2 0 -F[5 ;1 ] ffi [7]/=3 0 F[g], (20) 

where F indicates the number of times the correspond¬ 
ing state (I) appears in the product. Moreover, the sub¬ 
script of F represents the Young tableau for the state (I), 
■F[3,3]=5, -F[4,2]=9, E [5)1 ]=5, and F [6] =l. Consequently, 
there are 6 different kind of | J, S) states of the dibaryon 
which dictates the possible two baryon decay mode in 
strong force. In order to investigate the stability of the 
dibaryon against a strong decay, we will examine energies 
in relation to the threshold for the decay mode. 


TABLE III: The decay mode for the dibaryon into two 
baryons with respect to (I,S) states. 



AA 

AN 

AN 

NN 

AA 

NN 

(I.S) 

(3,0) 

(2,1) 

(1.2) 

(1,0) 

(0,3) 

(0,1) 


B. Color basis functions 

The dibaryon of our interest with isospin symmetry has 
color singlet function represented by the Young tableau 
of [2,2,2], as we mentioned earlier. Since the dimension 
of the Young tableau of [2,2,2] is five, there are five color 
singlet functions corresponding to the Young tableau of 
[2,2,2]. We define the color singlet functions as the fol¬ 
lowings : 



1 

2 

|Ci>= 

3 

4 


5 

6 


{[(12) 6 3] 8 [4(56) 6 ] 8 }i = 


(W (W (3 )q k (5 )e lmn q l (2 )q m (4 )q n (6) 
-^e ijk q i (l)q^2)q k (5)e lmn q l (3)q m (A)q n (6) 

- \e ljk q\lW (3 )q k (4 )e lmn q l (2 )q m (5 )q n (6) 

+ \e ijk q\lW(2)q k (4)e lmn q l (3)q m (5)q n (.6) 
-le ijk q\l)qi(2)q k (3)e lmn q\4)q m (5)q n (6)), 



1 

3 

1 c 2 )= 

2 

4 


5 

6 


{[(12)g3] 8 [4(56) 6 ] 8 } 1 = 


^ M (1)9 J (2 )q k ( 5)ei mn q l (3 )q m (4 )q" (6) 
-^e ijk q i (l)q^2)q k (4)e lmn q\3)q m (5)q n (6) 

+ ^e ijk q i (l)q^2)q k (3)e lmn q l (A)q m (5)q n (6)), 



1 

2 

|C 3 )= 

3 

5 


4 

6 


{[(12) 6 3] 8 [4(56)3] 8 } 1 = 


(W (W (3 )q k (4 )e lm nq l (2 )q m (%" (6) 
-\e ijk q i (l)q^2)q k (4)e lmn q l (3)q m (5)q n (6) 

+ ±e ijk q i (l)qi(2)q k (3)e lmn q\4)q m (5)q n (6)\ 


It is very important to understand the property of basis 
functions of the dibaryon in calculating the expectation 
values of both the confinement and hyperfine potential, 
proportional to either or A^A^CTi • <jj, respectively. 
So we will now establish the basis functions of the color, 
isospin, and spin based on the Young tableau, which are 
very useful for constructing completely antisymmetric 
states. Then, the expectation values can be easily cal¬ 
culated by using the complete antisymmetry properties. 



1 

3 

|C 4 )= 

2 

5 


4 

6 


{[(12) 3 3] 8 [4(56) 3 ] 8 }i = 


^=( e ijk q\l)qi(2)q k (4)e lmn q l (3)q m (5)q n (6 ) 
-y ijk q i (l)q^2)q k (3)e lmn q l (A)q m (5)q n (6)), 
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1 

4 

\Cs)= 

2 

5 


3 

6 


{[(12)g3] a. [ 4 ( 56 ) 3 ] i}i = 


^e ijk q\lW(2)q k (3)e lmn q l (4)q m (5)q n (6). 


( 21 ) 


We note that the color singlet functions followed by the 
standard Young-Yamanouchi representation is symmet¬ 
ric with respect to any adjacent particles that he in the 
same row, and is antisymmetric with respect to any par¬ 
ticles that lie in the same column. The definition next to 
the Young tableau expresses the convenient intermediate 
states for constructing the color singlet, and originates 
from the [8]<E>[8], and the [l](g)[l]. The orthogonality of 
the color singlet functions, (Ci\Cj)=5ij is easily obtained 
by using the tensor form, and in fact, results from the 
orthogonality of the standard Young-Yamanouchi bases. 

There are several ways for calculating the expectation 
value of XfXj with respect to the color singlet functions. 
Among those, it is very useful to consider one that is 
based on the irreducible matrix representation of the 
permutation group with respect to the standard Young- 
Yamanouchi bases whose irreducible matrix for the trans¬ 
position is symmetric. Moreover, when considering the 
operator for describing three gluon exchange, through ei¬ 
ther i fabcX^ Xj Xj , or dabcX^XjXj, as we shall show the 
detailed calculation in Appendix, this method gives us a 
simple form if we know the irreducible matrix represen¬ 
tation of the permutation group. However, in our case of 
the fully antisymmetric color <S> isospin £g) spin state, la¬ 
belled by | CiljSk), we only need the expectation value of 
AJAlj because the expectation value of A°A| as to a possi¬ 
ble as the expectation value of A^A^ can be obtained from 
the former using the antisymmetric property of the ba¬ 
sis functions. In fact, this calculation is performed using 
the formula, A°Aj=-8/3N, where N is the number 
of the participant particle in this dibaryon, N=6, result¬ 
ing in 15(AJA2)=-16. Here, 15 is the number of ways of 
pairing between particle i and particle j (i < j, i, j =1, 
2, 3, 4, 5, and 6 ). 


C. Flavor basis functions 

Since the flavor of the dibaryon is in the irreducible 
representation of SU(2), we consider the flavor basis func¬ 
tions in terms of the isospin representation that is allowed 
to the dibaryon. As in the case of color basis functions, 
the isospin part can be obtained using similar techniques 
based on Young tableau. Moreover, it is convenient to es¬ 
tablish the orthogonal basis functions with a certain sym¬ 
metry, making use of the standard Young-Yamanouchi 
bases. For the dibaryon, as mentioned above, the repre¬ 
sentation of isospin has 1=0, whose Young tableau is [3,3] 
with dimension of 5, and 1=1, whose Young tableau is 
[4,2] with dimension of 9, and 1=2, whose Young tableau 
is [5,1] with dimension of 5, and 1=3, whose Young 
tableau is [6] with dimension of 1 : 


• 1=0 : 

I I°i)= 

14 °)= 
• 1=1 : 
\ii)= 

1 1 \)= 

14 )= 
• 1=2 : 
\ll)= 


1 m- 


I4 2 )= 

• 1=3 : 

|/ 3 >= 


5 basis functions with Young tableau [3,3] 

14 °)= 


1 

2 

3 

4 

5 

6 


1 

2 

4 

3 

5 

6 


14 °)= 


1 

3 

4 

2 

5 

6 


1 

2 

5 

3 

4 

6 


14 °)= 


1 

3 

5 

2 

4 

6 


9 basis functions with Young tableau [4,2] 

14 )= 


1 

2 

3 

4 

5 

[ad 



1 

2 

3 

5 

4 

6 



1 

3 

4 

5 

2 

6 



1 

3 

4 

6 

2 

5 



14 )= 
I il)= 


1 

2 

3 

6 

4 

5 



1 

2 

5 

6 

3 

4 



14 )= 
14 )= 
14 )= 


1 

2 

4 

5 

3 

6 



1 

2 

4 

6 

3 

5 



1 

3 

5 

LelJ 

2 

4 



5 basis functions with Young tableau [5,1] 

1 - 4 ) = 


1 

2 

3 

5 


4 


1 

3 

4 

5 

[cd 


1 - 4 )= 


1 

2 

3 

4 

6 

5 


1 

2 

4 

5 

Leh 


1 basis function with Young tableau [6] 


1 

2 

3 

4 

5 

6 


We can also see the orthogonality of the isospin states, 
(/(|/j) = Sjj in a given irreducible representation of 
isospin from the orthogonality of the standard Young- 
Yamanouchi bases as well as the orthogonality between 
any two different irreducible representations, according 
to the group theory. 


D. Spin basis functions 


For the spin states of the dibaryon, the representa¬ 
tion of spin states of dibaryon is the same as that of the 
isosopin states because the irreducible representation of 
SU(2) should also be applied in this case : 


• S=0 : 

|S?) = 

|S 4 °) = 


5 basis functions with Young tableau [3,3] 


1 

2 

3 

4 

5 

6 


1 

2 

5 

3 

4 

6 


|S 2 °)= 

|sg)= 


1 

2 

4 

3 

5 

6 


1 

3 

5 

2 

4 

6 


IS 3 °)= 


1 

3 

4 

2 

5 

LelJ 


S=1 : 9 basis functions with Young tableau [4,2] 

I Sl)= 


1 

2 

3 

4 

5 

6 



|si)= 


1 

3 

4 

5 

2 

6 



l^)= 

l^)= 


1 

2 

3 

5 

4 

6 



1 

2 

3 

6 

4 

5 



|si)= 

\sk)= 


1 

2 

4 

5 

3 

6 



1 

2 

4 

6 

3 

5 




















































































































































































































































7 


|5]) = 


1 

3 

4 

6 

2 

5 



15])= 


1 

2 

5 

6 

3 

4 



|5])= 


1 

3 

5 

6 

2 

4 



• S=2 : 

I sf)= 
\S$)= 

|5f)= 


5 basis functions with Young tableau [5,1] 


1 

2 

3 

4 

5 

6 


1 

2 

3 

5 

6 

4 


1 

3 

4 

5 

6 

2 



|5f)= 


|5 4 2 )= 


1 

2 

3 

4 

L®lJ 

5 


1 

2 

4 

5 



• S=3 : 

|5 3 )= 


1 basis function with Young tableau [ 6 ] 


1 

2 

3 

4 

5 

L®lJ 


We are now in a position to construct the completely 
antisymmetric function of the color ® isospin ® spin 
state of dibaryon, which we will denote by |Cj, Ij, £*,). 
In particular, we choose the |/, £) basis function to be 
in the [3,3] representation, and in the conjugate of the 
color singlet [ 2 , 2 , 2 ], so that the color singlet <g) isospin 
® spins becomes fully antisymmetric. With this IS cou¬ 
pling scheme, we can find the fully antisymmetric func¬ 
tion of | C,I,S) for all of (I,S), by using the Clebsch- 
Gordan (CG) coefficient for making the representation 
of [3,3] of isospin <g) spin function. In calculating the 
Clebsch-Gordan coefficients of the element of the permu¬ 
tation group, Sq, it is convenient to use the factorization 
property that factorizes the CG coefficients of S n into an 
isoscalar factor, which is called K matrix, and a CG co¬ 
efficients of S „ i . For example, the isoscalar factor can 
be defined by [ 2 f|, 


sQfWv , [f'VW\[f\pov) = 

where S in the left-hand (right-hand) side is a CG coef¬ 
ficients of S n ( S n -i ). In this notation, [f p ] is the Young 
tableau associated to £ n _1 which can be obtained from 
[/], the Young tableau of S n , by removing the n-th par¬ 
ticle characterized by pqy, where p is the position of the 
n-th particle in the row, q , the position of the (n-l)-th 
particle in the row, y 1 the position of the (n- 2 )-th parti¬ 
cle in the row, respectively. In our case, by repeating the 
process of factorizing the CG coefficients of Sq further, 
we find the CG coefficients from the following formula, 


5 ( [f'jp'q'y'r' [f"]p"q"y"r" | [f]pqyr) = 
^([/V[/>"I[/]p)^([/;^ , [^]<z"|[/p]9) 

(23) 


where S in the third row is the CG coefficient of £ 3 . 
When we calculate the CG coefficients, we use the rel¬ 
evant isoscalar factors for £ 4 , £ 5 , and £r in Eq. (23) 
which is obtained by Stancu and Pepin [30] • In the case 
of (I,S)=(0,1), the | [Y°, £ x ]) basis functions belonging to 


the Young tableau of [3,3] are presented as the followings 


|[/°,5 1 ] 1 )= 

Vg 


1 

2 

3 

4 

5 

6 


9 

V 10 


A 0 ) ® l-Si) 


^ lr0 


9 


I A) ® IS 2 1 ) 


it!) ® IS 3 ) -g-|^ 2 ) ® |S!) 


v “K)®ls!> 


18 ; l©®isi)-f' 




Veo 


2 V 5 

— 

Vrn 

Vgo 


A 0 ) ® |£s)- 
A°)®|5])- 
AVI5])- 


a/30 , 


1 a 0 ) ® v + Via 0 ) ® v - Via 0 ) ® 15 ])+ 


lg I 4/ — I 36 

IAW 


18 


V60 0 \ ^ | cl 


36 


(24) 


IA°,£ 1 ] 2 )= 


1 

2 

4 

3 

5 

6 


vAO 1 r0\ 1 ci\ t/20 q 1 r0 \ I c i\ 


■ia°) ® 15 ]) - Via 0 ) ® 15 ]) 


18 

A 61 r o\ ^ 1 ci\ 


18 

+ ^1® ® is!) - ^ 


18 

A 


IA) ® |5g) 


IA°)®|5 2 1 ) + ^-|A°)®|5l) 


V20 lr0 , nfnU , vAo, r0N _ ,„ lv , V30 


18 


■IA) ® |5]) 


9 AW + 


+ ^ IA°) ® |5i) + ^IA°) ® 15]) + ^plA°)« |5g) 

- ^IA°) ® |5 2 ) - ^IA°) ® 15]) + ^IA°) ® |5]) 

+ ^IA°> ® |5]) + ^flA°) ® |5l> - ^|A°) ® 15]). 

(25) 


|[ j0 ,£ 1 ] 3 )= 


1 

3 

4 

2 

5 

6 


^IA°)®|5 4 V^IA°)®|5]) 


^, r o 


18 

AA 7 - 0 \ o I cl\ 


18 

vAo 


18 


A) ® |5g) 




IA°) ® |5]) X^IA 0 ) ® |5g) 


36 


V 6 lr0 . VW IT o^,„u vAO, 


+ ^IA 0 )®|5])- 18 


IA) ® |5]) - Vl J °) ® |5]) 


18 


t/lO g : vAO g 


~VI A)® 15])- g 


-IA) ® |5]) - Vl J °) ® |5|) 


36 


+ AA 0 ) ® |S 3 > - 4 P|iJ> ® l«r) - +Pl 4 °) ® IS!) 


18 

v/l5 

HT 


18 


IA°) ® |5 4 ) 


36 

' /I5, O®|s!) + ^|/ 5 0 )®|s!). 

(26) 


18 





































































































































\[I°,S%)= 
v 7 ^ 


1 

2 

5 

3 

4 

6 


18 

VIE 

HT 

VlE 

IF 

v 7 ^ 

~V2~ 


I4)®l4) + 
14°)® I 4) + 
I4°)®l^> + 


VTE it0v _ ,„iv v 7 ^ 


18 

-v/15 

HT 

v 7 ^ 


14°) ®I4V 

|/ 2 °)®|5 5 1 ) + 


18 

VEO 

FV 

V6, 


I4)®I4) 
14°)® I 4 1 ) 


36 l4°)®l^)-^-l4 0 }®l^) 


44)®I^ 1 )-^I4°}®|5 9 1 )- 


12 


(27) 


I [J°,5 1 ] b )= 

v 7 ^ 


1 

3 

5 

2 

4 

6 


-14°) ® 14) + 

+ 4^l4)®l4)- 


V15 it0x _ io1v x/15 


18 

v 7 ^, 


-|4)®I4)- 


+ 


36 

VlE 

IF 

F 

6 


18 

v/ ^Q I 7^0\ ICl\ 
18 


-I4)®I4) + 


l4)®l4>-4Pl4°)®l4)- 


-^I4°)®|5 1 1 )-^I4°)®I^ 1 )- 


36 
v 7 ^, 


18 

VIE 
1 F 

VEo 

FT 


14°) ®I4) 
14°) ®I4) 
I4)®l4) 
(28) 


Coupling the isospin Cg) spin basis function obtained from 
the IS scheme to the color singlet basis function, we find 
the color ® isospin ® spin state satisfying the fully anti¬ 
symmetry property for (I,S)=(0,1). This is given by 

14 1°, S 1 ) = -^(|Ci) ® \[I°,S%) - I C%) ® |[T°,5' 1 ] 4 >- 

|C3>®|[/ 0 ^ 1 ]3) + |C 4 )®|[/°,5 1 ] 2 )- 
ICs)®!)/ 0 ^ 1 ]!)). (29) 


In the case of (I,S)=(1,0), the If/ 1 , 4]) basis functions 
belonging to the Young tableau of [3,3] are presented as 
the followings : 


|[4,S°] 2 >= 


1 

2 

4 

3 

5 

L®J 


^14) ® I 4°) - ^14) ® 14°) + ^pl4) ® 14°) 
^14) ® 14°) - ^14) ® 14°) + ^14) ® 14°) 
—14) ® 14°) + ^14) ® 14°) + ^14) ® 14 


18 


^I4)®I4°) + ^P|4)®I4°) 
^14) ® 14°) + ^14) ® 14°) + ^14) ® I-4) 
^ 14) ® I4°) - ^pi4) ® 14°) + ^14) ® 14°) 
^14) ® 14°) + ^14) ® 14°) - ^pl4) ® 14°)- 

(31) 


|[4,4] 3 >= 


1 

3 

4 

2 

5 

6 


- ^14) ® I-S'?) - ^14) ® 14°) + ^fl4) ® 14°) 

- ^14) ® 14°) - ^14) ® I -S'?) - ^14) ® I-S'?) 
+ ^14) ® I s’?) - ^14) ® l-s 3 °) - ^pl4) ® I-S'?) 

- ^14) ® 14°) - ^14) ® 14°) - ^fl4) ® 14°) 

\/l5 


■ 


^14) ® 14°) - ^Pl4) ® 14°) - ^14) ® 14°) 
^14) ® 14°) + 4) ® 14°) + ^14) ® 14°)- 

(32) 


14,4]!) = 



1 

2 

3 



4 

5 

6 



^I4)®l4)- 
^14)® 14°) 


^I4)®I4°) + 


18 13 

4 io 


g-|4) ® 14°) + 


18 

v 7 ^ 

~F 


14) ® 14°) - ^ 14) ® 14°) + 
14) ® 14°) + ^Pl4) ® 14°) - 


^I4)®I4 0 )- 

^I4)®I4 0 )- 

^14) ® 14V 

VEo 


36 


18 


^Pl4)®l4 0 )- 


I4)®I4°)+ 

(30) 


I4,5°] 4 )= 


1 

2 

5 

3 

4 

6 


-#4>®I4V 
-^I4)®I4°) + 
-^|4)®I4°) + 
+ ^fl4)®l4 0 )- 


^I4)®I4°)-^?I4)®I4°) 

VIE 


18 


18 

v 7 ^ 

"36" 


14) ® 14°) + ^?I4) ® 14°) 


36 

46, 

6 


14) ® 14°) - V4) ® 14°) 


® |s 6 °>- 


12 


(33) 
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I4,s°] s )= 


1 

3 

5 

2 

4 

6 


^Vi> ® is?) 

V30, r i 




18 

a/30 


+ -^14) ® i^ 2 °) - Vi 4) ® is 3 °> + Vi4) ® \s° 3 ) 


18 

4l5._ 1 . 


36 
Vl5| ri . 


18 

430, rl . 


+ -^414) ® 1*4) - Vl4) ® 1*4) - Vl4) ® 14) 


18 
y/30 . 1 . 


18 


36 


12 


- ^1+ ® is?) - 14') ® |Sg>- 


(34) 


Likewise, we find the color 0 isospin 0 spin state satis¬ 
fying the fully antisymmetry property for (I,S)=(1,0) to 
be given by 


I C,I\S°) = -^(|C7i) ® \[I\S°] 5 ) - | C 2 ) 0 |[/\S°] 4 )- 

\C 3 )0\[I\S O } 3 ) + \C 4 :)0\[I 1 ,S O } 2 )~ 
|c 5 )® |[ J T 1 ,5'°]i>). (35) 


In the case of (I,S)=(1,2), the jf/ 1 , S' 2 ]) basis functions 
belonging to the Young tableau of [3,3] are presented as 
the followings : 


I4,S 2 ]i) = 

y/TE 


1 

2 

3 

4 

5 

6 


4i5. 


14) ® 1*4) + Vl J 2> ® 14) + 44) ® |S 2 >- 


^14) ® I 4) - \\ll) 0 |S 5 2 ) + ^\ll) 0 15?)+ 


^i4)®i4) 

42 


^14) ® 14)-^14)® 14)- 


g 1 ^ 7 ) ® |<S’|) g~l4) ® 1*4) — l-fg) ® |S|). (36) 


|[/ 1 ,s 2 ] 2 )= 

y/15 


1 

2 

4 

3 

5 

6 


I 4) ® 14) - d4) ® 14) + Vl4) ® 14)- 


\/l5| r i. 


^14) ® 1*4) + ^1 4) ® 1*4) - ^14) ® |S|>- 


^1^,0+ , V ° , rl \ <e»Jlo2\ , VX " U -|J 6 1 )® |S|)- 


9 H4)®l4) + ^l4)®l4) + ^°' 


""g“l4) 0 |Sg) + -|/g) ® |S 2 ) — g 14 ) ® I'S' 2 ) - 
^14) ® 14) - ^14) ® 14) + ^14) ® |s 5 2 ). (37) 


I 4,S 2 ] 3 )= 


1 

3 

4 

2 

5 

6 


^14) ® 14) - ^14) ® 14) - ^14) ® 14)+ 
^14) ® 14) - ^14) ® 14) - ^14) ® 14)- 
^ 14) ® 14) - §14) ® 14) + ^14) ® 14)+ 
^14) ® 14) - ^14) ® 14) - §14) ® 14)+ 
^I4>«14) - ^14) ® 14) + ^14) ® 14)- (38) 


|[4,s 2 ] 4 >= 


1 

2 

5 

3 

4 



^pl4) ® I 4) + ^14) ® 14) + ^14) ® 14)- 
^14) ® 14) - ^14) ® 14) - ^14) ® 14)— 
^14) • 14) + ^14) ® 14) + ^14)•I4>- 


4270, rl 


45 


14)® 14)- 


(39) 


|[4,s 2 ] 5 )= 


1 

3 

5 

2 

4 

6 


4® I rl\ „ I rr2\ I 243, r i^| C 2\ 


^14) ® IS 2 ) - 2^14) <8 |4) + -g 


-14) ® 14)- 


46, r i 

9 


4e lr i v 

9 


14) ® 14) - ^rl4) ® I 4) + ^14) ® 14)- 


^14) ® 14) + ^14) ® 14) + ^14) ® 14)- 


4270, 1 


45 


14)® 14)- 


(40) 


We find the color 0 isospin 0 spin state satisfying the 
fully antisymmetry property for (I,S)=(1,2) to be given 
by, 


144, S 2 ) = -J=(|Ci) ® |[4,S 2 ] 5 ) - |C 2 ) ® |[7 1 ,5' 2 ] 4 >- 

|C' 3 )®|[/ 1 ,S 2 ] 3 ) + |C 4 )®|[/ 1 ,S 2 ] 2 )- 
ICs^IIAS 2 ]!)). (41) 


In the case of (I,S)=(2,1), the |[7 2 , S 1 ]) basis functions 
belonging to the Young tableau of [3,3] are presented as 
the followings : 
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l[/ 2 ^ 1 ]l}= 

s/15 


1 

2 

3 

4 

5 

6 


x/l5, 


I ID ® \sl) + VI J 2 2 ) ® l^ 1 ) + d4 2 ) ® &)- 


\\il) ® l^ 1 ) - lui) ® |5i) + ^|/i 2 ) ® |5 5 1 )+ 

^1 It) ® ISb 1 ) + wpl 7 !) ® I*-*]) - ^I4 2 ) ® |S]>- 
^I4 2 ) ® 1$) - II) 0 15s 1 } - ^|/f) 0 |S]). (42) 


|[/ 2 ^ 1 ] 2 )= 
a/15 


1 

2 

4 

3 

5 

6 


VT5, 


|/|) ® itf) - dif) ® l^ 1 ) + V|/ 2 2 ) ® 1^}- 


^1) ® l^ 1 ) + ^ I il) ® 1$) - ^|/f> ® is])- 

^ \ll) ® I*]) + ^|i?) 0 I'S'g) + ^|J|) 0 |<S])— 
^I4 2 ) 0 l<S]) + \\ll) 0 |S]) - ||Z|) g |$)- 

^1 II) 0 ISl) - ^|/1) 0 15s 1 ) + ^|/f) 0 |S]). (43) 


1 

3 

4 

2 

5 

6 


IMs)-^ 

^I4 2 ) 0 Itf) - i|/|) 0 15a 1 ) - ^|/f) 0 |S 3 1 )+ 

^pl il) 0 isl) - \\ll) 0 \Sl) - ^\ll) 0 |S])- 

^I4 2 ) 0 I-?]) - ^|/ 5 2 ) 01 'Stf) + ^1 il) 0 \s$)+ 

^fl ii) 0 1 $) - ^I 4 2 ) 01-57 1 } - \\il) 0l$)+ 

^I4 2 ) ® l^ 1 ) - ^|/ 3 2 ) ® l^g 1 ) + ^|/ 4 2 ) 0 l^g)- (44) 


|[ 4 2 , 5 1 ] 4 )= 


1 

2 

5 

3 

4 

6 


0 \Sl) + ^|/ 3 2 ) ® \Sl) + II) 0 |S£>- 
^1 il) 0 IS]) ^ ^|if> 0 IS]) - ^|/ 3 2 ) 0 |S])- 

\/3, r2v n 1 , a/ 3 , r2x „ 1 ^,1 x , v^5 


9 -1^)0 IS]) 
Iff 


9 I4 2 )0|s 7 1 ) + ^|/i 2 )0|s 8 1 )- 


-|/|) 0 |s]). 


(45) 


|[Z 2 ,S 1 ] 5 )= 


1 

3 

5 

2 

4 



^1© ® I© - ^1 © ® I© + ^pl© ® |SJ>- 
^1© ® I© - ^ © ® I© + ^ I © ® I©- 
^1© ® I© + ^1© ® is?) + ^1© ® is?)- 


/270 

45 


I II) 0 IS]). 


(46) 


We find the color ® isospin ® spin state satisfying the 
fully antisymmetry property for (I,S)=(2,1) to be given 
by, 


\C, 1 2 , S 1 ) = -^(|Ci) 0 |[/ 2 , S X ] B ) - |C 2 ) 0 |[/ 2 , S 1 ],)- 

|C' 3 )0|[/ 2 ,S 1 ] 3 ) + |C 4 )0|[/ 2 ,5 1 ] 2 )- 
|C 5 )0|[/ 2 ,S 1 ] 1 )). (47) 

For the case of (I,S)=(3,0) and (I,S)=(0,3), we find 
straightforwardly the fully antisymmetric color ® isospin 
0 spin state, written by, respectively, 

|C,/ 3 ,S°) = -^(ICi) 0 |/ 3 ) 0 |S°) - |C 2 ) 0 |/ 3 ) 0 |S2)- 

v 5 

| C 3 ) 0 |/ 3 ) 0 |S 3 °) + |<7 4 ) 0 1I 3 ) 0 |S 2 °)- 
|C' 5 )0|/ 3 )0|5 1 0 )), (48) 

\C,I°,S 3 ) = -^(|Ci) 0 |/ 5 °) 0 |S 3 ) - |C 2 ) 0 |/°) 0 |S 3 )- 

|C 3 ) 0 14°) 0 IS 3 ) + IC 4 ) 0 |I 2 °) 0 |S 3 )- 
IC 5 ) 0 |lf) 0 |S 3 )). (49) 

In dealing with the expectation value of —A^A j(Ti • aj 
with respect to | C,I,S) state for all of (I,S), the sym¬ 
metry property of the state make this calculation sim¬ 
ple in that ' ct : ,)=15(-A);A 2 (Ji • er 2 ) as ar¬ 

gued early. Moreover, using the symmetry properties of 
the wave function, one can derive the effective formula 
which is expressed in terms of the Casimir operators of 
the isospin, spin, and color, given by |31], 

-E^A CA 7 CfT * ■ = 

[|jV(JV - 6) + 4/(7 + 1) + p(s + 1) + 2cy, (50) 

where N is the total number of quarks in this system, and 
C c =iA c A c , that is, the first kind of Casimir operator of 
SU(3) in the system of N quarks. Since we must consider 
only the color singlet as a physical observable, the term 
C c vanishes. For practical purposes, our the calculation 
of (—A 3 A 2 <7 i • er 2 ) can be easily performed with the sym¬ 
metry of between particles 1 and 2 in the |C, I, S) state, 
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whose the property of the symmetry is definitely derived 
from the Young tableau. Then we find the following for¬ 
mula: 


| -> = A?A§ | 


CTl • (72 


>=- 


). CTi • <7 2 



)• (51) 


TABLE IV: The expectation value of —A^cr; • aj with 
respect to \C,I,S) for all possible (I,S) configurations in the 
dibaryon, which will be denoted by Vd■ AV is ^-(VIi+V^), 
in which the Vbi and Vb 2 are those of the baryons to which 
the dibaryon can decay. 


(I.S) 

(3,0) 

(2,1) 

(1.2) 

(1,0) 

(0,3) 

(0,1) 

v d 

48 

80 

3 

16 

8 

16 

8 

3 

AV 

32 

80 

3 

16 

24 

0 

56 

3 


V. NUMERICAL RESULTS 


TABLE V: The mass of the dibaryon in (I,S) state with the 
two types of potentials given in Eq. (2) and Eq. (3). The 
binding energy Eb is taken to be the difference between the 
mass of the dibaryon and the two baryon threshold. The 
dimension of the variational parameters are given in fm -2 . 


(i,s) 

(3,0) 

(2,1) 

(1,2) 

(1,0) 

(0,3) 

(0,1) 

Type 1 

3132.3 

2926.5 

2808.5 

2710.8 

2808.5 

2639.6 

Variational 

o=1.2, 

a=0.9, 

a=1.4, 

a=3, 

a=1.4, 

a=3.6, 

parameters 

6=1.2, 

6=1.5, 

6=1.8, 

6=1.4, 

6=1.8, 

6=1.4, 


c=1.4 

c=1.7 

c=0.9 

c=1.3 

c=0.9 

c=1.3 

E b 

597.3 

681.9 

563.9 

756.6 

273.5 

685.4 

Type 2 

2845.5 

2711.7 

2629.9 

2558.6 

2629.9 

2504.1 

Variational 

o=0.6, 

a=0.5, 

a=0.7, 

0=2.2, 

a=0.7, 

a=2.6, 

parameters 

6=0.7, 

6=0.9, 

6=1.1, 

6=0.8, 

6=1.1, 

6=0.8, 


c=0.6 

OO 

o 

II 

o 

c=0.5 

c=0.7 

c=0.5 

c=0.8 

E b 

371.1 

498.2 

416.4 

606.0 

155.5 

551.5 


In this section, we analyze the numerical results ob¬ 
tained from the variational method, by using the com¬ 
pletely symmetric spatial function as the trial function. 
TablelVlshows the result of the analysis with the trial spa- 
cial wave function given in Eq. (ITSl) after adding 14 addi¬ 
tional forms as discussed before. Among all the dibaryons 
with (I,S), it is the dibaryon with (1=0,S=3) that is most 
likely to be stable against the strong decay. However, as 


we see in the Table 0 , we find that even for this state 
the two baryon threshold lies below the dibaryon mass. 

Although simple comparison of ' a j) 

the hyperfine interaction between the dibaryon and two 
baryons, as given in Table. ITVl shows the splitting is min¬ 
imal in the (1=0,S=3) channel, the lowest mass of the 
dibaryon for both types of the potential considered is far 
above the threshold of two A baryons. 

Table m shows that the dibaryon mass is larger for 
type 1 potential compared to that for type 2. This is 
because the size of the dibaryon in the former is smaller 
than that in the latter, as determined by the inverse of 
the variational parameters a, b, c. To better understand 
this point, we show the values of each energy term of the 
dibaryon with (1=0,S=3) in Tabic IVII 


TABLE VI: The values of each energy term of the dibaryon 
with (1=0,S=3) and A baryon. AE is the difference between 
the dibaryon and two A baryon in each term. 


Type 1 

kinetic 

linear 

coulomb 

hyperfine 

Dibaryon 

1282.8 

2603.6 

-446.9 

186.9 

A Baryon 

Variational parameters 

589.2 

1214.8 

-239.3 

111.6 

o=1.4 6=2.1 

AE 

104.2 

173.8 

31.6 

-36.3 

Type 2 

kinetic 

1/2 power 

coulomb 

hyperfine 

Dibaryon 

722.7 

2924.1 

-346.6 

132.5 

A Baryon 

Variational parameters 

333.2 

1410.4 

-186.6 

81.7 

o=1.4 6=0.7 

AE 

56.2 

103.2 

26.7 

-30.8 


Since the size of the dibaryon in both types of poten¬ 
tials are larger than that of a single baryon, the value of 
the kinetic part of the dibaryon is comparatively larger 
than the sum of that of the two A baryons. Moreover, 
the value of kinetic part of dibaryon in type 2 is a little 
larger compared to that in type 1 due to its relatively 
small increase in size. 

In addition to the kinetic term, the effect of 1/2-power 
confinement part of the type 2 potential causes the mass 
of the dibaryon to decrease compared to the case of type 
1. Nevertheless, the lowest mass in type 2 is still about 
155 MeV above the threshold of two A baryons, and no 
other choice for the confinement potential, such as 1/3- 
power, is expected to change the stability of the dibaryon. 

In investigating the stability in the present work, 
we can consider three-body color confinement operators 
which are mentioned in the introduction, as this approach 
may change the stability of multiquark configurations, 
such as that of the dibaryon. The two types of operators 
can be expressed in terms of the permutation operators 
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given by 

d a b c F a F b F c + (»fcj)] + 1/ 

- + (* fc ) + O'*0]> 


f<*°F?FjF£ = -±[{ijk)-{ikj)], (52) 

where I is the identity operator, and ( ijk ) and (*j) 
are operators belonging to the permutation group of 
Se, which is called 3-cycles, and 2-cycles, respectively, 
and F“=1/2A“. Also, there is another formula for 
d abc F^F b F^ : which can be conveniently shown to be in¬ 
variant under the SU(3) algebra, given by [27j . 


d abc F^ Fj 


g \yi+j+k 2 ^i+j+k 



(53) 


where C^ is the first kind of Casimir operator, and C^ 2 \ 
the second kind of Casimir operator of SU(3). Since 
baryon consists of three quarks, the SU(3) invariant op¬ 
erators are written by, 


d abc F?F b F£ =^[(123) + (132)] + 

-±[(12)+ (13)+ (23)], 


fabcpapbpc = _ * [( l 23 ) _ ( 13 2)]. (54) 

For baryon which has one color singlet represented by 
the standard Young-Yamanouchi basis of Young tableau 
[1,1,1], we use the irreducible representation with one 
dimension given by 



1 


1 

a 

2 

= (-ir 

2 


3 


3 


(55) 


where (— 1) G is 1 if a is an even permutation, and -1 if 
a is an odd permutation. By using this formula, we can 
easily calculate the action of the operators on the color 
singlet of baryon, \C)=^=tij k q l {l)q> {2)q k (‘i)) in Eq. (10) 


d abc F?F b F£\C) = y| C), (56) 

f abc F?F%F£\C) = 0. (57) 

As we can see, while the introduction of d ahc F{ 1 F. b F£ can 
contribute to the baryon mass, / abc i 7 ' 1 0 F|i 7, | will not. 

Likewise, for dibaryon which has five color singlet bases 
corresponding to the standard Young-Yamanuochi bases 
of Young tableau [2,2,2], we can calculate the irreducible 
matrix representation of the permutation operators be¬ 
longing to Sa in the form of 5 x 5 matrix in terms of the 


five color singlet bases. Since the irreducible matrix of 
(ij) and [(ijk) + (ikj)] are symmetric, and the irreducible 
matrix of (ijk) is antisymmetric, d abc F?F b F£ has a form 
of symmetric matrix and f abc F?'F b Ffc has a form of her- 
mitian matrix with vanishing diagonal elements. As we 
will show in detail in appendix, the matrix of d abc F'jF, b F-j 
and f abc F° F. b F-j are given by, 

(Cild^FfFlF^ICi) = (C 2 \d abc F?F b F£\C 2 ) = 

(C 3 1 d abc F® F$ F£ | C 3 ) = (C 4 \d abc F?F%F*\C 4 ) = ~ 

(C 5 \d abc F?F$F£\C 5 ) = 10/9, 

(Ci | f abc F® F% F 3 1 Cr) = (C 2 \f abc F?F b F£\C 2 ) = 
(C 3 \f abc F?F b F£\C 3 ) = (C 4 \f abc F?F b Fg\C A ) = 
(C 5 \f abc F?F b F§\C 5 )=0. (58) 

Then, we can find the expectation value of the three-body 
confinement operators, d abc F? F|F 3 , and f abc F{ l F, b F^ in 
terms of |C, I 1 , S '- 7 } for (I=i,S=j), reminding that 

\c,r,s j ) = ±={\Cx) 0 |[r,^] 5 ) - |C 2 ) 0 |[r,^] 4 )- 

|c 3 ) 0 |[/ i ,^] 3 ) + |c 4 ) 0 |[r,^] 2 )- 
|C' 5 >0|[/\^] 1 }), (59) 

one finds the followings: 

(C,r,Si\d abc F?F b FZ\C,P,Si) = ±(-4 x A + H) 


{C, P, S j \f abc F?F b F£\C, r, S j ) = 0. ( 61 ) 

Due to the complete symmetry of | C,I l , S J ) under 
any permutation of Sq, one finds that ( C,I\S J | 
d abc F t a F^F^ | C,I\S j ) =( C,I\Si | d abc F?F$F§ | 
C, P , S J ) for any l, m, n (/<m<n=1^6). Similar relation 
holds for ( C, I\S j | f abc Fj 1 F% l F£ | C, I\ S j ). As far as 
the color factors are concerned, one finds that the contri¬ 
bution of the d-type three body interaction contributing 
to the dibaryon as given by Cf = 20 times Eq. (I6U1) is 
equal to twice that of the baryon given in Eq. (15B1) . 

If we attempt to add the d-type three-body confine¬ 
ment potential, the result will depend on the detailed 
functional form multiplying the SU(3) invariant opera¬ 
tors. For the simplest choice, we can chose it to be the 
sum of two body potential such as ( Vi 23 = d abc F 1 Q E 2 b F| ( 
^ 12 /uo +i"i 3 /ao + ^ 23 /ao) )• However, for such a simpli¬ 
fied form, the result will only be a re-parametrization of 
the parameters of the model and there will be no change 
in the stability of the dibaryon. On the other hand, a 
intrinsic three body force will change the situation and 
on that grounds, it is of great importance to have some 
ideas on the understanding of the two-body and three- 
body confinement. 
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VI. SUMMARY 


as 


M 


In order for the total wave function to be fully anti¬ 
symmetric in six quarks system with only u, d quarks, 
we first consider the spatial function that is fully sym¬ 
metric, and find 15 Jacobian coordinates appropriate for 
the symmetry, and then, construct the spatial function 
desirable for this scheme with Gaussian spatial function 
to perform variational method in nonrelativitic Hamil¬ 
tonian. Secondly, we classify the physical states with 
respect to isospin (I) and spin (S), and find the color 
singlet basis functions, isospin basis functions, and spin 
basis functions allowed to the six quarks system, and con¬ 
struct the color ® isospin <g) spin states that should be 
completely antisymmetric, by means of IS scheme that 
couples the color basis function to IS basis function. We 
find that there does not exist a compact dibayron system 
in all system that is stable against the decay into two bay- 
rons with corresponding quantum numbers. Hence, the 
recently observed peak in the / = 0, S = 3 B=2 chan¬ 
nel [mi should be a molecular configuration composed 
of two A dibaryons. 


Appendix A 


(123) = (23) • (12) = (1/ + + \k A§) = 

1 1 3 

-I + - KXj + 2 d abc F?F%F% + 2if abc F?F%F£. 

i<j 

(A2) 


(132) = (23) • (13) = (±7 + ±A£A§)(!/ + TjA^Ag) = 

1 1 3 

-I + - ^ A) : A) + 2 d abc F?F%F£ - 2if abc F?F%F£. 

i<j 


(A3) 


By adding Eq. (A3) to Eq. (A2) and subtracting 
Eq. (A3) from Eq. (A2), we obtain d^F^F^F^ and 
j abc papbpC' Also, we can apply this process to any 
(ijk) (i<j<k=1^6), and finally obtain the Eq. (52) ; 


ja-bc zpa ppb 


rpa Tpo rpc 
“i “j “k 


= + {ikj )] + 

- ^[(u) + (*&) + (Jk)], 


In this section, we derive the three-body color opera¬ 
tors, which is invariant to SU(3), in terms of the relevant 
permutation operators. This will enable us to represent 
the the three-body color operator with respect to color 
singlet basis functions of the dibaryon. We can express 
the algebra of SU(3) as the permutation of two particles, 
as in the case for the algebra of SU(2) ; 

(1.2) = ^01 -cr 2 , 

(I2) = h+\\ c l \ c 2 - (Al) 

Here, I is the identity operator, and (12) is 2-cycles per¬ 
mutation. Then, noting that (123), and (132) can be 
written as (123) = (23) • (12) and (132) = (23) • (13), 
we can straightforwardly present 3-cycles permutation 

I 


f abc F?F b F c k = --[{ijk) - {ikj)]. (A4) 

Now, we can construct the matrix representation of 
the three-body color operators in terms of the standard 
Young-Yamanouchi bases corresponding to the color sin¬ 
glet bases of dibaryon. As mentioned earlier, the stan¬ 
dard Young-Yamanouchi bases which is orthonormal to 
each other are written by, 



(A5) 


The matrix of d abc F°-F b F k is given by, 


/ 

36 


d abc F?F b F£ = 


0 - 
0 
0 
0 


1 0 

— 0 

36 u 

0 — — 
u 36 


0 
0 
0 

0 36 

0 0 


0 \ 
0 
0 

0 
10 
9 / 


,dabcF?F$F£ = 


( -- 

36 

0 
0 
0 
0 


0 0 

-M 0 
0 


0 

0 

0 

35 

36 
5 

9V2 


o \ 

0 

0 

5 

9\/2 

0 J 
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(~ A 


dabcF?F%Fg = 


36 

0 


0 

0 

0 

0 

25 

0 

5 

5 

36 

6\/3 

6\/6 

0 

5 

36 

0 

0 

5 

0 

5 

5 

6\/3 

36 

18\/2 


V 0 6/6 


0 - 


18/2 


,dabcF?F$Fg = 


( - A 
36 

0 
0 
0 
0 


0 

0 

0 

0 

25 

0 

5 

5 

36 

6\/3 

6V6 

0 

5 

36 

0 

0 

5 

0 

5 

5 

6\/3 

36 

18V2 


V 0 eSff 


0 - 


18/2 


d ahc F*F*Fi = 


36 

0 

0 

0 

0 

\ 


( 

35 

72 

5 

8y/3 

5 

8y/3 

5 

24 

5 

12V2 

0 

5 

0 

0 

0 




5 

5 

5 

5 

5 

36 




8V3 

72 

24 

24V3 

12V^ 

0 

0 

25 

36 

5 

6\/3 

_ 5 
6y/6 


ri 1 F a F b pc — 

1 ^abc 1 l i 3 i 5 — 


5 

8\/3 

5 

24 

5 

72 

5 

24y/3 

5 

12\/6 

0 

0 

5 

5 

5 




5 

5 

5 

5 

5 

6V3 

36 

18V2 




24 

24 V3 

24V3 

72 

36\/2 

V 0 

0 

_ 5 
6\/6 

5 

18V2 

0 

) 


\ 

5 

12V2 

5 

12V6 

5 

12\/6 

5 

36\/2 

0 


/ 35 5 _5__5_ __ 


d abc F?F»Fg = 


72 

8\/3 

8^3 

24 

12V2 


72 

8^ 

8^ 

24 

12\/2 

5 

5 

5 

5 

5 


5 

5 

5 

5 

5 

8\/3 

72 

24 

24^/3 

12V6 


8\/3 

72 

24 

24^/3 

12V6 

_5_ 

8V3 

5 

24 

5 

72 

5 

24V3 

5 

12^ 

,d abc F?F*Fg = 

_5_ 

8y/3 

5 

24 

5 

72 

5 

24V3 

5 

12V6 

5 

5 

5 

5 

5 


5 

5 

5 

5 

5 

24 

24V3 

24V3 

72 

36\/2 


24 

24V3 

24 

72 

36\/2 

, 5 

\ 12V2 

5 

12V6 

5 

12V6 

5 

36\/2 

0 J 


, 5 

\ 12V2 

5 

12V6 

5 

12V6 

5 

36V2 

0 

/ 35 

72 

_5_ 

8V3 

5 

8\/3 

5 

24 

5 

12V2 

\ 

/ -- 

36 

0 

0 

0 

0 \ 


d ab cF?F*FC = 


_5_ 

8\/3 

5 

8\/3 

_5_ 

24 

5 


_5_ 

72 

_5_ 

24 

5 

24\/3 

5 


_5_ 

24 

_5_ 

72 

5 

24a/3 

5 


24\/3 

5 

24\/3 

_5_ 

72 

5 


12>/6 

5 

12>/6 

5 


12^2 12>/6 12\/6 36\/2 


36\/2 

0 


,da b cF?F»F* = 


_ _5_ 

36 

0 

0 

0 


0 

25 

36 

_5_ 

6/3 

5 

6/6 


0 

5 

6/3 

_5_ 

36 

5 

18/2 


0 

5 

6/6 

5 

18/2 

0 


d abc F^F^Fi = 


/ - A 
36 

0 
0 
0 
0 


_ _5_ 

36 

0 

0 

0 


0 

0 

0 

0 

25 

5 

36 

6V3 

5 

5 

6V3 

36 

5 

5 

6\/6 

18V2 


0 

0 

5 

6/6 

5 

18/2 

0 


) d a bcF£ F3F5 = 


t 35 
72 

_5_ 

8\/3 
5 

8\/3 

_5_ 

24 
5 


_5_ 

8\/3 

_5_ 

72 

_5_ 

24 
5 

24-\/3 

_ 5 

12>/2 12->/6 


5 

8\/3 

_5_ 

24 

_5_ 

72 

5 


24 V3 

5 _ 

12V6 36\/2 


_5_ 

24 

5 

24^3 

5 

24a/3 

_5_ 

72 
5 


12y/2 

5 

12\/6 

5 


12y / 6 

5 

36\/2 

0 


7 papbpc _ 

Ll abc r 2 r 3 r 6 — 


dabcF2F»Fg = 


35 

5 

5 

5 

72 


8-s/3 

8/3 

24 



5 

5 

5 

8V3 

72 

24 

24/3 



5 

5 

5 

8V3 

24 

72 

24/3 

5 


5 

5 

5 

24 


24 V3 

24/3 

72 

5 


5 

5 

5 

12V2 

121/6 

12/6 

36/2 

/ 

35 5 

5 

5 



72 8/3 

8/3 

24 



5 5 

5 

5 



8/3 72 

24 

24/3 



5 5 

5 

5 



8/3 24 

72 

24/3 



5 5 

5 

5 



24 24/3 

24/3 

72 



5 5 

5 

5 


V 

12/2 12/6 

12/6 

36/2 


12 V 2 

5 

12 V 6 

5 

12 V 6 

5 

36\/2 

0 


12V2 

5 


, dabcF£F%F£ = 


/ 35 

72 
5 

8\/3 

_5_ 

8V3 

_5_ 

24 

_3— 

\ 12V2 


,d abc F2F*FC = 


/ -A 

36 

0 

0 
0 
0 


_ _5_ 

36 

0 

0 

0 


5 

5 

5 


5 

8/3 

8/3 

24 


12/2 

5 

5 

5 


5 

72 

24 

24/3 


12/6 

5 

5 

5 


5 

24 

72 

24/3 


12/6 

5 

5 

5 


5 

24/3 

24/3 

72 


36/2 

5 

12/6 

5 

12/6 

5 

36/2 


0 

0 

0 

0 

\ 


0 

0 

0 



25 

5 

5 



36 

6/3 

6/6 


7 

5 

5 

5 



6/3 

36 

18/2 



5 

5 

0 

) 


6/6 

18/2 
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dabcF^F^ = 


36 

0 

0 

0 

0 

\ 


/ 

5 

36 

0 

0 

0 

0 

\ 

0 

25 

0 

5 

5 




0 

25 

0 

5 

5 


36 

6V3 

676 




36 

673 

676 


0 

0 

5 

36 

0 

0 


,dabcF£F*FZ = 


0 

0 

5 

36 

0 

0 


0 

5 

0 

5 

5 




0 

5 

0 

5 

5 


673 

36 

18 72 




6 73 

36 

1872 


V 0 

5 

6\/6 

0 

5 

18V2 

0 

) 


l 

0 

_5_ 

6\/6 

0 

5 

18V2 

0 

/ 


d abc Fq = 


( -A 

36 

0 

0 
0 
0 


0 

36 u 

0 - — 

u 36 


V 


0 
0 
0 

35 5 

36 9v/2 


0 \ 
0 
0 


5 

9V2 


0 


> d a bcF£ F§ Fq = 




/ -A 

36 

0 

0 
0 

V 0 


0 0 
-M 0 

0 -m 


o \ 
o 
o 


0 -M 0 


0 0 


12 
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(A6) 


We can prove that Y^i<j<k dabcFp F 1 - F£ is invariant to SU(3), which means that it commutes with the generators of 
dibaryon, F a = 1/2 ( A“ + Aj + Ag + Af + A!? + Ag ), as it is proportional to the identity operator as can be seen 
in the following: 


E dabcFf Fj F£ = 
i<j<k 


( f 0 0 0 0 \ 

0 f 0 0 0 

0 0 f 0 0 

0 0 0 f 0 

Vo o o of/ 


(A7) 


The matrix of fabcF^F’-F^ is given by, 


( o 


f Z7 a T?b Z? c _ 

Jabc-T 1 -C2 -^3 — 


n/3 

4 

0 

0 

0 


i7 3 
4 

0 

0 

0 

0 


0 

0 

0 

n/3 

4 


0 0 \ 

0 0 

i7 3 n 
4 u 

0 0 
0 0 / 


f JP a rpb jpC _ 

7 Jabc^l *2 ^4 — 


/ o 

0 

0 0 

0 0 - 

V o o 


0 0 
0 0 
0 


i* T? a T?b jpc _ 

Jabc-T i r 2 r 5 — 


V- 


0 

0 

0 

_ i_ 

4 

i 

2s/2 


0 0 
0 - 


4 

0 - 


— i — \ 

2 V 2 

0 


2%/3 


0 - 


t 0 
0 


2V6 


f jpa rpb rpc _ 

Jabc-Fl -^ 3-^4 — 


V 


0 

i7 3 
4 

0 

0 

0 


n/3 

4 

0 

0 

0 

0 


0 0 

0 0 

0 - 


2\/3 2\/6 

0 0 

0 0 

o \ 

0 


f J7'a Tpb rpc _ 

7 Jabc-F 1 ^2 -^6 — 


4V3 

i 

76 

0 


0 


4%/3 

i 

2 V 6 


4 ^ 

0 

2 

2V2 


2\/6 

i 

272 

0 


r z?° z?fr z? c _ 

, Jabc^l -C3-C5 — 


2\/2 

0 

0 

0 


0 - 

0 

0 


4 

0 - 


4y/3 

0 

0 


0 \ 

0 

_ i 

~ 76 

0 


0 / 

_ i_ i \ 

4 9./9 


4 

0 - 

i 

2^3 


4 2\/2 

0 0 

i i 

2V3 2V6 

0 0 


2 V 6 


_ i i7 3 
\ 4y/2 4 V2 


2^3 
i i 
4\/6 4^2 


0 / 


_ i_ i \ 

4 4'/2 

is/3 


0 - 


472 


2y/3 4Ve 

0 - 


4 V 2 

0 


r z? a z?& z? c _ 

J abc-T l -C 3 -Tg ~ 


t 0 
0 
0 

A 


0 

0 

i_ 

4 

0 - 

iV 3 _ 

4\/2 


0 


_i i \ 
4 4y/2 

i7 3 

4^2 


0 


273 

i i 
476 4V2 


2y/3 4y/E 
0 - 


4y/2 

0 


( 0 
0 


f ipa z pb rpc _ 

7 Jabc^l -^4-^5 — 


0 -i 
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(A8) 


In this case, it turns out that the SU(3) invariant operator is fabcFiF%F£ + fabcFiF^F^ + f a b c F^F^F§ FfabcFiF^FQ 

1 £ TP& T?b 77'C 1 £ rpa irb 77'c 1 i* ZT'a 7716 z^c 1 i 1 rpa irb 77'C 1 i* rpa rpb 77'C 1 x rpa rpb 77'C i* 771a rpb 77'c 

“T Jabc-F l -^3 -^4 + Jabc-F 1 -^3^5 + /a6c-^l-^3^6 + Jabc-^l *4 -^5 + Jabc^ l-^4^6 +/a6c-^l -^5 -^6 " Iabc ^2 ^3 P 4 ~ 

f , rpa zpb 77'C ^ , rpa T^b tt'c £ , 771a 77*6 77'C £ , ZTia 7716 77'C i 1 , rpa rpb tt<c £ , jr<a rpb t?c £ , rpa rpb rpc £ , rpa rpb 77'C 

Jabc? 2-^3-^5 " J abc^ 2 ^3 ^6 ~ Jabc-^2 ^4 ^5 ~ Jabc^ 2 ^ 4^6 ~ Jabc^ 2-^5-^6 " Jabc^ 3^4^ 5 " Jabc-F 3 ^4^6 " Jabc-F 

- fabcF£F§ Fq , due to the fact that this operator is 


/o 0 0 0 o\ 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
V 0 0 0 0 0 / 


(A9) 
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